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i ' The dominant source of decoherence for an electron spin in a quantum dot is the hyperfine interaction with the surrounding bath of 
' nuclear spins. The decoherence process may be slowed down by subjecting the electron spin to suitable sequences of external control 
^ ' pulses. We investigate the performance of a variety of dynamical decoupling protocols using exact numerical simulation. Emphasis is 
given to realistic pulse delays and the long-time limit, beyond the domain where available analytical approaches are guaranteed to work. 
Our results show that both deterministic and randomized protocols are capable to significantly prolong the electron coherence time, even 
when using control pulse separations substantially larger than what expected from the upper cutoff frequency of the coupling spectrum 
between the electron and the nuclear spins. In a realistic parameter range, the total width of such a coupling spectrum appears to be the 
physically relevant frequency scale affecting the overall quality of the decoupling. 



S. 

1 Introduction 

Q ■ The dynamics of electron spins in semiconductor structures, such as quantum dots (QDs), has attracted 
, CJ, . much attention due to both its fundamental importance, as well as the potential applications in spin- 
based quantum information processing: not only does the electron spin degrees of freedom provide a 
. natural two- level system with long energy relaxation time [1] , but semiconductor-based quantum computing 
[ architectures can benefit from well-developed technology and a high degree of scalability [2j. Practical 
I/-) ' exploitation of such a potential, however, is severely hindered by the loss of phase coherence caused by the 
unavoidable coupling between the electron spin and its surrounding environment. For typical QD operating 
£^ | conditions and geometry, decoherence mechanisms associated with phonon processes or magnetic spin-orbit 
■ interactions are effectively suppressed, leaving the hyperfine coupling with nuclear spins as the dominant 
, decoherence channel. The resulting decoherence dynamics occurs very rapidly, with the coherence times 
from ensemble measurements (TJ ) being of the order of tens of nanoseconds for typical GaAs QDs at 
experimentally relevant (sub-Kelvin) temperatures and moderate (sub-Tesla) magnetic fields [3]. 
As a result, viable methods for decoherence suppression are of key relevance. A number of proposals have 
' been put forward for increasing the coherence time, for instance by seeking a high degree of nuclear spin 
C , polarization 0H], by narrowing the nuclear spin distribution (6j[7], or by cleverly disentangling the bath 
from the electron spin via external control [8]. However, these methods remain challenging to implement 
by current or near-term capabilities, and alternative approaches are still actively sought. 

In this context, a promising approach is offered by dynamical decoupling (DD) methods [9]. Borrowing 
^ ' their inspiration from coherent averaging techniques in high-resolution NMR spectroscopy |10|lllj. DD 
^ . schemes aim at removing the effect of environmental coupling by subjecting the electron spin to suitably 
tailored sequences of control pulses. The long correlation time and distinctively non-Markovian behavior 
of the nuclear spin reservoir make the DD approach naturally suited for electron spin coherence control 
in QDs, as demonstrated by the recent implementation of single-pulse Hahn echo experiments [3|ll2j. 
which succeeded at extending Tjf by two orders of magnitude. Thus, it is both timely and important 
to undertake a comprehensive theoretical study of performance of multi-pulse DD protocols, in order to 
assess relative strengths and weaknesses, and identify optimal schemes. Only a few protocols have been 
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quantitatively analyzed so far. In addition, the majority of investigations have focused on the case of QDs 
subjected to strong magnetic fields, predicting a coherence enhancement of at least an order of magnitude 
in GaAs QDs [13] . Still, the experimentally important situation of low-to-moderate fields is significantly 
more difficult to analyze due to the competition between purely adiabatic (energy-conserving) decoherence 
and relaxation processes. Existing studies have mostly relied on restrictive approximations - in particular, 
semi-classical treatments [14] or Magnus expansion (ME) techniques [15]. Thus, little is known about the 
actual performance of DD protocols in a variety of practically important regimes. 

Our goal in this paper is to further push a quantitative analysis of the electron spin decoherence problem 
in regimes which are relevant for experimental DD implementations in QDs, yet have received marginal 
consideration so far [16J. Special attention is devoted to the limit of zero external field, where dephasing 
and relaxation must be simultaneously eliminated, and the asymptotic long-time limit, where the effects 
of finite control rate and coherent error accumulation are crucial, so that neither the ME nor the standard 
quasi-static approximation for the nuclear spin bath are reliable a priori. 

The content is organized as follows. Section [2] provides the essential background on the QD model Hamil- 
tonian and various DD protocols. Section [3] summarizes both the relevant steps of the methodology used 
to tackle the problem, and presents a comparative analysis of all the DD protocols under investigation. 
A detailed discussion of the best DD protocol based on periodic concatenated design is presented, includ- 
ing the effects of control pulse delay, level of concatenation, intra-bath dipolar interaction, and coupling 
spectrum. Concluding remarks are given in Section [4] 

2 Quantum dot model and dynamical decoupling protocols 

The dynamics of an electron spin S coupled to a bath B of N nuclear spins in a QD is described by a 
total Hamiltonian of the form 



Here, S and denote the electron and the fc-th bath spin operators, respectively; ojq is the electron 
Zeeman energy in an external magnetic field Bq; H$b gives the hyperfine contact interaction between 
the electron spin and the nuclei; and Hb corresponds to the intra-bath dipolar coupling between nuclear 
spins |17U18][T9lll2] . The parameters and T^i determine the strength of the respective interactions. In 
particular, the coupling energy per nucleus depends on the electron probability density at the nucleus 
location. We focus on the limit of zero external field Bq = 0, so that all directions are equally affected by 
decoherence, and choose = 1/2, since the analysis does not depend qualitatively on the specific spin 
value (up to an appropriate renormalization of the parameters). The Zeeman energies of the nuclear spins 
and the anisotropy of the hyperfine coupling are very small and can be neglected for the time scales we 
are interested in. The bath is assumed to be initially unpolarized, which is justified by the fact that in 
typical experiments the temperature is much larger than the characteristic energy of nuclear spins. 

In the absence of control pulses, the electron spin described by the above Hamiltonian undergoes free 
induction decay (FID) on a time scale given by T 2 * = (iV^ 2 /8)- 1 /2 ; where A = (T, k A l/ N ) 1/2 

is approx- 
imately KTVeV for typical GaAs QDs with N = 10 6 [L9\. In our analysis, we set h = 1, and measure 
frequency and energy in units of A, with time correspondingly being given in units of 1 /A. 

In order to slow down the decoherence process, we recur to DD methods. The basic idea consists of 
modifying the system dynamics by subjecting the electron spin to a sequence of 7r^ pulses, where n denotes 
the control axis. In the simplest setting, the pulses are equally separated by an interval r, during which the 
system evolves freely. We assume bang-bang controls, that is, arbitrarily strong and instantaneous pulses, 
which we draw from a control group Q = {gj}, j = 0, . . . , \Q\ — 1, according to Pjj = gjg\. For a single spin 
as relevant here, the sequences are derived from the irreducible Pauli group Qp = {I , X,Y, Z} [9], where 
/, X, Y, Z are the identity and the Pauli matrices of the electron spin, respectively. A variety of control 
protocols have been developed, which exploit deterministic [HHIH] and stochastic [20l[2Tj design. In cyclic 



H = H,s + Hsb + Hb 
Hs = luqSz, Hsb = S • J2k=i ^fclfc, Hb = 
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DD methods, the pulse sequence is fixed, pre-determined, and periodically repeated in time; in randomized 
DD, the future control actions are not known in advance, but chosen at random from a set of alternatives. 

The construction of deterministic DD sequences has been mainly based on the average Hamiltonian 
theory [11], having as a goal the cancellation of the dominant terms in the ME series for the underlying 
time evolution operator. The simplest periodic DD (PDD) sequence we consider here ensures the removal 

of the zeroth order term TC of the average Hamiltonian TC at every T n = nT c = n\Q\r, where n £ N, T c 
is the cycle time, and \G\ = 4 is the size of the DD group. By fixing the path of Qp as {I,X,Y,Z}, the 
control sequence becomes C\ = CqXCqZCqXCqZ , where Cq denotes the interval of free evolution. Thus, 
the cycle propagator for a system subjected to the PDD sequence is given by 

U(T C ) = Texp j-i J °[H + H c (s)]ds\ = exp [-iHT c ] = exp -i(TC {0) +TL (1) + • • -)T C (3) 

= ZUqXUqZUqXUq , 

where T indicates time ordering, H c (t) corresponds to the applied control Hamiltonian, Uq = exp(— iHr) 
is the free propagator, and the zeroth order term in the ME is 

n {0) = Y,gtUo(t k+1 ,t k )g k = - [ZHZ + YHY + XHX + H] = H B , (4) 
' ' k 

that is, terms coupling the system and bath degrees of freedom are removed to lowest order. The terms 

depending on Hgp appear in the higher-order corrections, in particular TC ~ Ti oc t, which guarantees 
DD between the system and the bath in the limit of arbitrarily fast control, r — > 0. Clearly, r is necessarily 
finite in realistic settings, its minimum value being constrained by practical limitations (e.g., for QDs and 
other devices operating at dilution refrigerator temperatures, the need to avoid excessive radio-frequency 
heating of the sample, which imposes limitation on the total absorbed power). In order to improve DD 
performance, going beyond what is achievable by the first-order schemes above, two main strategies are 
available: higher-order DD, ensuring cancellation of higher order terms in the effective Hamiltonian; and 
reduction of error accumulation in time due to the remaining terms of the Hamiltonian. 

As prototypes of higher-level deterministic sequences, we consider symmetric DD (SDD) and con- 
catenated DD (CDD) |15j . The SDD sequence corresponds to a symmetrization of C\, being writ- 
ten as Cs = CqXCqZCqXCqICqXCqZCqXCq. Its cycle time is twice as long as PDD, but at every 

T n = 2nT c = 8nr, besides Ti^° , all odd terms in TC are also canceled, leading to TC oc r 2 . CDD relies on a 
temporal recursive structure, so that at level £+1 the pulse sequence is given by Ci + i = CgXCgZCiXCiZ . 
Here, we truncate the concatenation procedure at a certain level and repeat a periodic sequence, referred 
to as PCDD^, after every 4^r. For example, I = 2 leads to PCDD2, which is also a symmetric protocol, 
although usually superior than SDD at reducing error accumulation. 

The main goal of stochastic protocols is to guarantee that errors accumulate probabilistically in time, 
instead of coherently as in deterministic sequences. As representatives of these methods, we consider naive 
random DD (NRD), random path DD (RPD), and symmetrized random path DD (SRPD). NRD corre- 
sponds to a sequence of pulses chosen randomly from Qp with equal probability. For RPD, randomization 
is associated with the path to traverse the group Gp, that is, at every T n = 4r one chooses at random one 
of the following options: {C XC YC XC Q Y , CqXCqZCqXCqZ, CqYCqXCqYCqX , CqYCqZCqYCqZ , 
CqZCqXCqZCqX , CqZCqYCqZCqY}. SRPD improves over RPD by symmetrizing the random paths as 
in SDD and inserting a random pulse between paths. Besides slowing down the error accumulation at 
long times, RPD and SRPD also ensure good performance at short times, since they lead to effective 
Hamiltonians H e fj oc r and H e jf oc r 2 , respectively. 

Note that for randomized DD the periodicity, as well as the notion of a control cycle, is lost in general. 
Thus, the comparison between deterministic and randomized protocols is most naturally carried out in 
the so-called logical frame. In this frame, H c {t) is removed from the evolution through the transformation 
U(T) = Uc(T)U(T), where tilde indicates logical frame, U C (T) = Texp[— i J Q T H c (s)ds] corresponds to the 
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control propagator, and U(T) is the propagator in the physical frame, as given by Eq. ([3]). For deterministic 
sequences, both physical and logical propagators stroboscopically coincide at the completion of each cycle, 
but the same is not necessarily true for the stochastic methods. In the latter case, in order to return the 
state of the system to the physical frame, it suffices to keep track of the applied control pulses. 



3 Performance of control protocols 
3.1 Control metric 

To compare the above-mentioned control protocols in a way which avoids dependence on the initial state 
of the systems, |^s(0)), we shall consider the worst-case scenario and compute the minimum input-output 
pure-state fidelity, 

F m (T) = mm\ M0)) Tr\ps(T)pg(p)], (5) 

where ps{T) corresponds to the reduced density matrix of the electron spin at time T. Equivalently, F m (T) 
may be thought of as the minimum state overlap between the intended and the decoupled system evolution. 

Analytical bounds on the worst-case pure-state expected fidelity have been established for various DD 
protocols [20j[2TJ[l5] , however they apply only to short evolution times and small pulse separations. Nu- 
merical simulations are mandatory to investigate DD performance at long times and for a broad range 
of r values. We are especially interested in pulse delays beyond the restricted domain where the ME 
is guaranteed to converge, oo c T c <C 1, where oj c denotes a high frequency cutoff. While in most situa- 
tions lo c is determined by Hb, for the QD problem the nuclear spin bath is to a very high accuracy 
non- dynamical, thus oj c is determined by the electron-nuclear coupling Hamiltonian H$b- This leads to 
uj c ~ Ylk l-^fcl/4 ~ NA/4, which is very large for typical QDs. We will consider r ~ l/2er, where a repre- 
sents the power spectrum width of the entire system, with 2a ~ (Y2k ^l) 1 ^ 2 = v^/V^ [22] . Thus, the basic 
cycle time T c = 4t ~ y/Nuj^ 1 will be up to a factor y/N longer than the formal (worst-case) convergence 
requirement of the ME. We shall also assume that > are uniformly distributed random numbers. 

The minimum fidelity is numerically obtained by employing the quantum process tomography 
method [23J as follows. The initial state of the entire system is taken as a direct product of the electron spin 
state |?/>s(0)) and the bath spin state \iPb(0)). Four different |^s(0)) are considered: | |), | |), (| |) + | |))/V2, 
and (| T) + i\ l))/V2] while a single \^b(0)} is assumed, corresponding to a random superposition 

IV'-b(O)) = z2i=i c%\4>i) of all possible tensor products of the form \<p) = | T)i ® | 472 ® \ J73 <8 • • • ® \ T)jv, 
where Cj are uniform random numbers. The time-dependent Schrodinger equation for the joint system plus 
bath dynamics is then solved by applying the Chebyshev polynomial expansion method to the evolution 
operator [18]. At final time T, a partial trace over the bath is performed, and the four resulting reduced 
density matrices are used to compute the 4x4 x- m &trix |23j . The ^-matrix is a superoperator that de- 
scribes the effective dynamics of the electron spin and allows for the reduced density matrix of an arbitrary 
electron initial state to be evaluated from the equation 

3 

PS{T)= K mPs (0)Kl Xmn (T), (6) 

m,n=0 

where Kq = I, K\ = X, K2 = —iY, and K% = Z. Since an arbitrary initial state of S in the Bloch sphere 
may be expressed as |^ s (0)) = cos(0/2)| |) +sin(0/2)e^| |), 9 £ [0,vr], tp £ [0,2tt], F m (T) simplifies to 

F m {T) = min^) (^(0)|p 5 (T)|^(0)) . (7) 

In practice, after determining the matrix Xi we obtain F m (T) from numerical minimization, by comparing 
the results over a statistically meaningful set of initial guesses for (9,<p), and then selecting the worst case. 
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3.2 Results 

Fig. dHa) shows F m (T) for the above-mentioned DD protocols. For comparison, we plot also the FID 
signal, for which no control is applied. Because the characteristic time scale td for the nuclear dipolar 
dynamics (resulting from Hb) is at least two orders of magnitude slower than the one due to H$b in 
typical QDs, setting Hb = is justified over the time scales studied here. As seen, all schemes lead to 
substantial enhancement of the electron spin coherence, provided that the control time scale is sufficiently 
fast compared to an effective upper cutoff determined by the spectral width 2a. Similar to the behavior 
observed for DD oil/ f noise [24J , this indicates how general convergence arguments may place too stringent 
a bound for realistic systems, actual performance depending sensitively on spectral details. 

Among the various protocols, NRD shows (not surprisingly) the worst performance, consistent with the 
expectation that its advantages over deterministic DD appear only when the control group is very large. 
PDD and SDD are, respectively, outperformed by RPD and SRPD, confirming the reduction of error 
accumulation associated with randomization. Interestingly, the RPD protocol provides worse fidelity than 
SDD at times smaller than 10, but eventually also outperforms SDD at longer times. The best performance 
for this system, however, is achieved by the deterministic PCDD2. Compared with the case of closed systems 
under the action of reducible control groups investigated in [21], this confirms how concatenated design 
is especially efficient for irreducible averaging. A feature which also proves beneficial for concatenation is 
the absence of internal bath dynamics |15j . 

100, 

90 

E 80 

^70 
60 
50 

10" 1 10° 10 1 10 2 10 3 10 1 10 2 10 3 
Time Time 

Figure 1. Minimum fidelity vs. time in the logical frame. The Hamiltonian parameters are Bo = 0, To = 0, and TV = f5. (a) 
Comparison of various protocols with r = 0.1. For NRD and FID, data points are acquired after every t, for deterministic DD (PDD, 

SDD, PCDD2), right at the completion of each cycle, while for RPD and SRPD this is done after every 4r and 8r, respectively. 
Random protocols are averaged over 10 2 control realizations, (b) Comparison of PCDD2 (solid lines) and PCDD4 (crosses) for t = 0.3, 

0.4, 0.5 from top to bottom. 

The above findings prompt a more in-depth analysis of the PCDD protocol. In Fig. [T](b) , we compare 
the performance of two levels of concatenation, £ = 2,4, for different values of r. As expected, the quality 
deteriorates as r increases but, unexpectedly, PCDD4 underperforms PCDD2 for larger r, suggesting that 
predictions based on applications of the ME [15] may not hold in this case. 

The dependence of the PCDD2 performance on the pulse delay r is investigated in Fig. 2(a), by focusing 
on the time T go % where the protocol attains -^(290%) = 0.9 for given r. Two main features are observed: 
(i) ^90% increases monotonically as r decreases - consistent with the fact that no crossover from decoher- 
ence suppression to decoherence acceleration is expected as long as r is sufficiently small with respect to 
the inverse spectral width of the environmental coupling 0123]. That is, in terms of the formal analogy 
established between decoherence control methods based on DD pulses and frequent quantum measure- 
ments [25 3 26] , only the counterpart of quantum Zeno effect (QZE) is observed in the relevant parameter 
regime, with no occurrence of "anti-Zeno" behavior; (ii) The PCCD2 performance depends quite sensi- 
tively on r, moderate changes in the value of r (say, by a factor of two) changing the coherence time T 90 % 

(2) 

by up to three orders of magnitude. For PCDD2, the leading non-zero term in the ME is 7i , so that 
the error amplitude per cycle is bounded by 0(|£/| 3 (<tt) 3 ), and therefore the total error amplitude after a 
number T 90 %/(\G\t) of cycles grows as O(T9 %o" 3 (|^|r) 2 ). As a result, within the regime of validity of the 
ME, the time T 90 % at which the total error probability 1 — F m (T 90 %) reaches 10% is expected to increase 




6 



Wenxian Zhang, V. V. Dobrovitski, Lea F. Santos, Lorenza Viola, and B. TV. Harmon 




Figure 2. (a) T 90 % vs. t: Comparison of PCDD2 for bath sizes TV = 15 (squares), 17 (circles), and 25 (triangles); Hb = 0. The time 
Tgo% depends strongly on t, varying by three orders of magnitude when r is varied by a factor of two. (b) The effect of dipolar 
interaction between nuclear spins on PCDD2 performance, F m vs. time for r = 0.25; (c) Tg % vs. r. For both (b) and (c), the 
parameters are TV = 15, Bo = 0, and uniformly distributed random numbers r^; £ [—To, To]; To = (solid line), To = 0.1 (dashed line). 

as 1/r 2 with decreasing r. The steeper dependence shown by the data in Fig. 2(a) indicates that a naive 
extrapolation of the above estimate in the long-time regime of relevance is not directly viable. 

Finally, we investigate the effect of interaction between the bath spins, taking the intra-bath coupling 
parameters in (pQ) as uniformly distributed random numbers in [— ro,ro], and assuming that the bath 
spins are located on a 3 x 5 piece of a two-dimension square lattice. One may expect that the effect of 
intra-bath coupling becomes important once the characteristic coupling strength Tq is comparable with 
the system-bath coupling A. However, Figs. EJb) and (c) show that the effect of intra-bath interactions 
becomes clearly visible already for much smaller values of Tq = 0.1 A, leading to a decrease of the coherence 
time Tgn.% by an order of magnitude in comparison with the case To = 0. Although the fast bath dynamics 
regime is not directly relevant to standard GaAs QDs, where Tq ~ 0.01^4, our observations may be useful 
as a first step towards characterization of DD in other systems. 

3.3 Coupling spectrum 

Additional physical insight into the performance of DD schemes may be gained by considering the coupling 
spectrum of the QD Hamiltonian (JTJ) in the absence and in the presence of control. A typical spectrum for 
uniformly distributed random Ak (obtained numerically) is shown in Fig. EJa). For the homogeneous case 
of constant A^ = A (below, we set A = 1 for convenience), the spectrum may be evaluated analytically, the 
result being also shown in Fig.[3|a). Physically, the limit A^ = A is equivalent to a semiclassical treatment 
of the bath, where the quantum- mechanical operator B = vlfclk, representing the Overhauser field 
exerted by the bath on the electron spin, is approximated by a randomly oriented classical vector B with 
zero mean and variance a [14 t ll9 ^ [T2] . 

For A^ = A, the operator of the total bath spin I = 1^ is an integral of motion, and the Hamiltonian 
(HJ) acquires a simple form of two coupled spins, 

H = AS • I = H diag + H off = A SJ Z + A(S X I X + S y I y ) , (8) 

where I 3> 1. Let us consider a state of the whole system | f; I, I z ), which corresponds to the electron spin 
up, and the bath having total spin / with z-projection equal to I z . Since the value of I z + S z is also an 
integral of motion, this state is coupled only to one state, | J.; I, I z + 1). Thus, the total Hamiltonian matrix 
can be separated into a number of uncoupled 2x2 blocks, each block involving a pair of levels | f; /, I z ) 
and | |; J, I z + 1). Correspondingly, the whole system can be viewed as a collection of uncoupled two-level 
systems, with state space spanned by basis states | f; /, I z ) and | J.; I, I z + 1). The energy splitting between 
these basis states, e = I z + 1/2, is determined by the diagonal part of the Hamiltonian ([8]), H^ag, whereas 
the coupling matrix element between the basis states is determined by the non-diagonal part H ff. 
The coupling spectrum G(u>) quantifies the density of bath states weighted by the corresponding coupling 
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Figure 3. (a) Spectrum of the QD for Bq = 0, To = 0, and N = 10. The bars are obtained for a single realization of random A^, 
whereas the solid line corresponds to the homogeneous limit, = A. (b) Pictorial representation of the coupling spectrum G(lj) (solid 

line) and the control pulse spectrum (dashed lines). 



strength. The spectral density n(io) of two-level systems having energy splitting e = u may be computed 
by a combinatorial argument: this is the density of the bath states having a given value of I z (since e = I z 
for I z 3> 1), which leads to 

n(w) = -jL- exp(-co 2 /2a 2 ). (9) 
V27T a 



The modulus square of the matrix element V which couples the basis states | j; I, I z ) and | I; I, I z + 1) is 
\V\ 2 = (A 2 /4)(I 2 — I 2 ), and it should be averaged over all two-level systems having the same I z (i.e., the 
same to) but different values of I. We then find that \V\ 2 = a 2 /2, and the coupling spectrum G{oj) for the 
electron spin plus bath is (to a good accuracy) a sum over all coupling spectra of the individual two-level 
systems [27] . 

G(f>) = \V\ 2 n(co) = -^exp(-a; 2 /2^). (10) 



Control modifies, in general, the above coupling spectrum. Under DD, in particular, the relevant spec- 
trum becomes a convolution of the original coupling spectrum G(u) with the spectrum of the control 
pulses (Fig. Eth)), which is given by a series of 5- like peaks separated by a distance of the order of 1/r [25J. 
The two central peaks of the controlled spectrum have the weight of 4/-7T 2 [25], while the contribution from 
the remaining peaks is less than 20% and thus may be omitted. As a result, in the presence of DD pulses, 
the decay is caused predominantly by the part of G{oS) which has u ~ tt/t. Such a contribution decreases 
exponentially, as exp(— uj 2 /2a 2 ), with decreasing r. Correspondingly, DD becomes efficient for r ~ 1/cr, 
the inverse width of the spectrum, not the inverse spectrum cutoff frequency io c (which is larger than a 
by a factor of y/N 3> 1). 

The fast decrease of G(u) with u> may provide an intuition as to why T 90 % depends so drastically on r, as 
seen in Fig. 2(a). If we assume, in analogy with [24 , 25J , that the time scale T go % is inversely proportional 
to the contribution of G(u>) at u> = qir/r (where q is some protocol-dependent parameter of the order of 
unity), then a strong dependence on r is indeed expected: 

T 90% ~^exp(gV/2crV). (11) 

Also, Eq. (|10p further demonstrates why the anti-Zeno-like behavior (acceleration of decoherence by DD 
pulses) is not observed in our simulations: G(lo) monotonically decreases with uj, and even for large T c the 
total power of the modified coupling spectrum cannot increase. Since, as noted, our results for N = 15, 
17, and 25 (Fig. 2(a)) are very close to each other, we expect that the basic qualitative features of our 
analysis remain applicable to realistic situations with large N. 
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4 Conclusion 



We have characterized quantitatively the performance of various DD protocols for an electron spin coupled 
to a nuclear spin bath, including both standard cyclic and randomized DD design. Two main points emerge 
from our analysis: on one hand, DD methods demonstrate a clear potential for significantly enhancing the 
electron spin coherence time for an arbitrary initial state, actual performance depending on both control 
and physical parameters. In particular, the frequency spectrum of the environmental coupling proves 
fundamental for identifying the minimum time scale requirements that must be met in order for DD 
to be effective. On the other hand, care must be taken in extrapolating simplified approaches (such as 
the ME) to parameter regimes outside the established domain of applicability, in particular for coherence 
preservation over long times. While additional investigation is needed to gain a complete picture (including, 
for instance, the effect of realistic control imperfections), it is our hope that the steady progress witnessed 
in experimental solid-state control techniques [28 , 29J may soon allow to validate the usefulness of DD in 
realistic QD devices. 
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